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We describe a new class of self-similar symmetric a-stable pro- 
■ cesses with stationary increments arising as a large time scale limit in 

a situation where many users are earning random rewards or incur- 
ring random costs. The resulting models are different from the ones 
studied earlier both in their memory properties and smoothness of 
the sample paths. 



CM 

& • 1. Introduction. With the dramatic increase of importance of commu- 

nication networks came the need to better understand their behavior at 
different scales. This requires a construction of stochastic models that can 
plausibly arise as the result of activities associated with such networks. Lim- 
| iting stochastic processes often scale, and one hopes that such models can 

■ provide insight into the scaling of properties of the networks. 

Perhaps the best known result of this type is the paper [21], where the 
limiting model turned out to be (depending on the relationship between the 
number of users and the time scale) either fractional Brownian motion or 
a Levy a-stable motion (this paper followed up and was an improvement 
of the earlier papers [38] and [36]). The fact that either a light-tailed but 
' ■^ | long-range dependent model or a heavy-tailed but short-range dependent 

model could appear has become an article of faith; see, for example, [7] 
for an application in a network context. Other heavy-tailed limiting models 
have appeared (see, e.g., [25]), but the limiting processes are not long-range 
dependent (more about this will appear in the sequel). 

In this paper we exhibit a natural situation where the limiting model 



belongs to a new class of a-stable models. It is a self-similar process with 
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stationary increments, and we will argue that the increments are long-range 
dependent. Let (Wk,k £Z) be a sequence of i.i.d. random symmetric vari- 
ables satisfying 

(1.1) F w (x):=P(W >x)~a$ v x- a 

as x — ► oo, where < a < 2 and aw > 0. Further, let (Vi, V2, ■ ■ ■) be a se- 
quence of i.i.d. mean zero and unit variance integer- valued random variables, 
independent of k E Z), defining a random walk S n = V\ + ■ ■ ■ + V n for 
n > 1. If one views 5 n as describing the "position" of the "state" of a user 
at time n, and VFfc the "reward" earned by, "cost" incurred by or "amount 
of work" produced by the user in state k, then the total reward earned by 
time n is 

n 

(1.2) R(n) = J2W Sj . 

j'=i 

Assuming that there are many such users earning independent rewards or 
generating independent work, it turns out that a properly normalized se- 
quence of rewards converges weakly to a limit, which we will call an FBM- 
1/2-local time fractional symmetric a-stable motion. This is a particular 
case of a larger class of models, FBM- //-local time fractional symmetric 
a-stable motion, where < H < 1 (these are self-similar with exponent of 
self-similarity H' = 1 — H+H/a). We will represent this process as a stochas- 
tic integral with respect to an a-stable random measure, with the integrand 
being the local time process of a fractional Brownian motion with exponent 
H (hence the name of the model). The increments of this process are gener- 
ated by a conservative null flow (see below for the details) and, hence, this 
process turns out to be different from all other classes of a-stable self-similar 
processes with stationary increments that have been considered so far in the 
extensive literature on the subject. 

Two remarks have to be made at this point. First of all, the only reason 
for assuming symmetry of (Wk,k £ Z) is that dealing with symmetric te- 
stable (SaS) models leads to simpler expressions and unified exposition for 
all < a < 2. Classes of nonsymmetric stable models parallel to those we 
are working with in this paper can be defined without difficulty, the case 
a = 1 being the exception. Under suitable tail conditions, the random reward 
scheme with appropriate translation and scaling will converge to these stable 
processes. Second, our processes are related to a family of limiting models 
obtained in similar circumstances (but with a single user) in [15]. In their 
case, the limiting process is self-similar with stationary increments, but not 
stable. 

This paper is organized as follows. In the next section we will summarize 
the required information on a-stable processes and random measures, on 
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self-similarity and on local times of fractional Brownian motions. Our process 
is formally introduced in Section 3. The properties of the increment process 
are discussed in Section 4. In Section 5 we study the smoothness of the 
sample paths of local time fractional stable motions through their Holder 
continuity properties. It turns out that local time fractional stable motions 
can be naturally written as sums of absolutely continuous self-similar stable 
processes, and the decomposition goes through the chaos expansion of the 
local times of fractional Brownian motions. This is done in Section 6. In 
Section 7 we prove the aforementioned convergence of the random reward 
scheme to the FBM-l/2-local time fractional stable motion. We conclude 
with some comments and a discussion of possible extensions in Section 8. 

2. Preliminaries. Throughout this paper, we will deal with SqS pro- 
cesses given in the form 

(2.1) X(t)= [ f(t,x)M{dx), teT, 

JE 

where T is a parameter space, M is a symmetric a-stable random measure on 
a measurable space (E,£) with a u-finite control measure m, and f(t,-) £ 
L a (m,£) for all t € T. See Chapter 3 of [34] for information on a-stable 
random measures and integrals with respect to these measures. 

If the parameter space is countable (e.g., if T = Z), and the process is 
stationary (under the usual left shift operator), then it has an integral rep- 
resentation as above, but the kernels /„, = /(n, •), n £ Z, are of a special 
form. Specifically, one can choose 

(2.2) /n(;E ) =an ( x )^__JL( x )j /o^(x), x£E, 

for n = 0, 1, 2, . . . , where (\>:E — > E is a measurable nonsingular map (i.e., 
a one-to-one map with both (f> and _1 measurable, mapping the control 
measure m into an equivalent measure), where 

0-n{x) = JJ UO (j^(x), X&E, 
j=0 

for n = 0,1,2,..., with u:E — ► {— 1, 1} a measurable function and /€ 
L a (m,£); see [29]. Many properties of the resulting stable process are closely 
connected with the ergodic-theoretic properties of the flow (the group of 
maps) (<j) n ,n £ Z), an important classification of which is into dissipative, 
conservative null and positive flows; see [29, 32] and [33]. In particular, a 
key idea in the latter two papers is that it is possible to view stationary sta- 
ble processes corresponding to dissipative flows as short memory processes, 
those corresponding to positive flows as infinite memory processes and those 
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corresponding to conservative null flows as processes with a finite but long 
memory. Good general references on ergodic theory are [1] and [17] . 

A stochastic process (Y(t),t > 0) is called self-similar with exponent of 
self- similarity H if, for all c > 0, the processes (Y(ct),t > 0) and (c H Y (t), t > 
0) have the same finite-dimensional distributions. Most commonly studied 
are self-similar processes with stationary increments [(Y(t + a) — Y(a),t > 0) 
has the same finite-dimensional distributions for all a > 0]. The common 
abbreviation for such a process is SSSI (self-similar stationary increments), 
or ff-SSSI, if the exponent of self-similarity H is to be emphasized. 

For SSSI processes with a finite mean, the exponent of self-similarity is 
restricted to the range < H < 1 (apart from degenerate cases) and, in that 
range, there is a unique if-SSSI Gaussian process. It has zero mean and 
covariance function 

Cov(Y(s),Y(t)) = EXM [t 2H + S 2H _ {t _ 8) 2H ]f 

< s < t. This process is called the fractional Brownian motion (FBM). 

In the a-stable case, < a < 2, the family of SSSI processes is much 
larger. The feasible range of pairs (a, H) is 

( , (0<H<l/a, if0<a<l, 

[ ' \0<H<1, ifl<a<2, 

and, apart from the case < a < 1 and H = 1/a, a feasible pair (a,H) does 
not determine the law of an SaS -ff-SSSI process. 

Remark 2.1. The class of SaS SSSI processes constructed in this paper 
has exponent of self-similarity in the range 

{\<H<l/a, if0<a<l, 
H = l, ifa = l, 

l/a<H<l, ifl<a<2. 

It has been a long-standing challenge to describe classes of symmetric 1- 
stable SSSI processes with H = 1 other than linear combinations of in- 
dependent symmetric 1-stable Levy motion and the straight line process 
Y(t) = tY(l),t > 0. The model developed in this paper provides, in the par- 
ticular case a = 1, an entire family of such processes. 

Two of the most well-known families of SaS -f/-SSSI processes (with < 
H < 1) are obtained by taking two of the many possible integral representa- 
tions of the fractional Brownian motion and modifying them appropriately 
(in particular, replacing the Brownian motion as an integrator by an SaS 
Levy motion). These are the linear fractional stable motion and the real 
harmonizable fractional stable motion. The linear fractional stable motion 
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belongs to the class of self-similar stable mixed moving averages described 
by [26, 27]. Its increment process is generated by a dissipative flow. On the 
other hand, the increment process of the real harmonizable fractional stable 
motion is generated by a positive flow. We refer the reader to Chapter 7 of 
[34] and to [10] for more information on self-similar processes. 

A fractional Brownian motion with any exponent of self-similarity < H < 1 
has a local time process (l(x,t),x G R,t > 0) that is jointly continuous in 
x and t [3]. The self-similarity property of the fractional Brownian motion 
immediately implies the scaling property of the local time process: for any 
c>0, 

(2.5) (l(c H x, ct),xeR,t>0) = (c 1 ~ H l(x, t), x £ R, t > 0), 
a somewhat more convenient form of which is 

(2.6) (h,(x,cb),x 6l,t>oj = ^Z^,tYxeR,t>oV 

It is a simple consequence of (2.5) and of Theorem 6, page 275 of [14] that, 
on a set of probability 1, 

(2.7) lim l(x, t)=oo for all xeR. 

t— >oo 

Similarly, the stationarity of increments property of the fractional Brow- 
nian motion implies a type of stationarity of the increments of the local 
time, which can be formulated as follows. Let (£1,3^, P) be the probability 
space on which the fractional Brownian motion and its local time process 
live. Then, abusing somewhat the term "law" by applying it to an infinite 
induced measure, 

the law of (l(x,t + h)(u) - l(x,h)(u),t> 0) 

(2.8) 

under P x Leb does not depend on h > 0. 

A modification of the proof of Theorem 1.2 that leads to Corollary 1.1 in 
[39] gives us that 

l(x,t) — l(x, s) 

(2.9 K:= sup -f-^ — ; ' ; NX „ <oo a.s. 

x£R (t-s) 1 - H (logl/(t-s)) H 

0<s<t<l/2 

and has finite moments of all orders. (Note that using instead the estimates 
in [6] gives a slightly worse power of the logarithm: H + 1 instead of H .) In 
particular, l(x,t) has moments of all orders finite and uniformly bounded in 
all real x and all t in a compact set. 
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3. FBM-iJ-local time fractional stable motions. We now introduce our 
class of models. Let (Q', !f,P') be a probability space supporting a fractional 
Brownian motion (Bjj(t),t > 0) with exponent of self-similarity H, and let 
/ = l(x,t) = l(x,t)(u') be its jointly continuous local time process. Let M be 
an SaS random measure on the space Sl'xK with control measure P' x Leb, 
where Leb is the Lebesgue measure on R. The random measure itself lives 
on some other probability space We define 



(3.1) 



Y(t) 



l(x,t)(u')M(duj',dx), 



t>0. 



Our first result below shows that (Y(t),t > 0) is a well-defined SaS process 
which is self-similar and has stationary increments. We call this process 
FBM-H -local time fractional symmetric a-stable motion. 

Theorem 3.1. The process (Y(t),t> 0) in (3.1) is a well-defined SaS 
process. It has stationary increments and is self-similar, with exponent 

(3.2) H' = 1- H + H/a = l + H( --1 

\a 

Proof. To show that Y is properly defined we need to check that 

/ r(x,t)(uj')-p'(du}')dx = E' [ l a (x,t)dx <oo. 
v ii Jr 



We have 



(3.3) 



E' / l a (x,t)dx= f E' 
Jr Jr 



l a (x,t)l( sup \B H (s)\ > \x\ 

0<s<t 



dx 



< [ (E'l 2 (x,t)) a/2 (p'( sup \B H (s)\ > \x\)) 1/9 dx 
Jr V \o<s<t J J 



with q = 1 — a/2. Since the moments of the local time are uniformly bounded 
and 

(3.4) / fp'fsup \B H (s)\ > \x\\) /9 dx<oo 

JR\ \0<s<t J J 

as the supremum of a bounded Gaussian process has Gaussian-like tails, we 
conclude that the left-hand side of (3.3) is finite and, hence, (Y(t),t > 0) 
in (3.1) is a well-defined SaS process. 

Notice that for any c > 0, k > 1, 6\, . . . , 6f. G R and ti, . . . ,tk > we have, 
using (2.6), 



( k 

Eexpl i^OjYictj 

V 7=1 



exp 



E' 



3=1 



d.r 
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exp 



E' 



i=i 



d.r 



«p(-«^->B-^g V (^, ti ) 

e X pf-c°( 1 -"> +H E' / 
\ Jm 

Eexp ^ OjC^^^Yitj)^ . 



Therefore, (Y(i),i > 0) is H '-self-similar, with H' given by (3.2). 

Furthermore, for any h > 0, k > 1, . . . , 0& G R and ti, . . . > we have, 
by (2.8) 

E exp ^ (Vfe + fr) - Y(h)) j 



(3.5) 



exp 



exp 



E' 



E' 



^20j(!(x,tj+h)-l(x,h)) 

3=1 



<7.r 



i=i 



dx 



Eex p^E% y (^))- 



Therefore, (Y(i),t > 0) has stationary increments. □ 
Remark 3.2. Observe that: 

1. For < a < 1 we obtain a family of if'-SSSI SaS processes with 
ff'€(l,l/a). 

2. For 1 < a < 2 we obtain a family of -ff'-SSSI SaS processes with 
ff'€(l/a,l). 

3. For a = 1 we obtain a family of 1-SSSI SaS processes. 

Notice that, for a^l, different choices of the fractional Brownian motion 
exponent of self-similarity H lead to a different exponent of self-similarity 
H' of the SaS process (Y(i),t > 0) and, hence, to a different process. On 
the other hand, for a = 1 the exponent of self-similarity H' is independent 
of i/. Nonetheless, the processes (Y(t),t > 0) corresponding to different H 
are different in this case as well, as will be seen in the sequel. 
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4. The increment process. An object of interest for an SSSI process is 
its increment process. It is a stationary process, and its memory proper- 
ties are often of interest. For example, the increment process of a fractional 
Brownian motion, the so-called fractional Gaussian noise, is a standard long 
memory (if H > 1/2) model that was used by Mandelbrot (see, e.g., [19, 20]) 
to explain the famous Hurst phenomenon. Similarly, the increments of the 
linear fractional stable motion are called linear fractional stable noise, and 
those of the real harmonizable fractional stable motion are called (real) har- 
monizable fractional stable noise. It is often believed that the properties of 
the fractional noises are largely determined by the exponent of self-similarity 
of the original process. One of the goals of this section (which studies the 
increment process of the FBM-Jf-local time fractional a-stable motion) is 
to shed some light on this question. 

Let, therefore, (Y(t),t > 0) be an FBM- if -local time fractional SaS mo- 
tion and consider its increment process 

(4.1) Z n = Y(n + l)-Y(n), n = 0,l,..., 

which will be called FBM-H -local time fractional SaS noise. 

A very important property of the FBM-7f-local time fractional SaS noise 
is given in the following result: 

Theorem 4.1. The FBM-H -local time fractional SaS noise is generated 
by a conservative null flow. 

Proof. Note that the FBM-if-local time fractional SaS noise has an 
integral representation 

(4.2) Z n = I { (Z(x,n + l)(u/) -l(x,n)(uj'))M(duj',dx), n>0. 

JW JR 

Let C be the space of continuous functions from R to R and a probability 
measure on C, under which the coordinate map is the fractional Brownian 
motion with exponent of self-similarity H. Let m be a u-finite measure 
on C defined by m = (P^ x Leb) o T , where T:CxK->C is given by 
T(lo',x) = ijj' — x, w'eC,x£l. Let L:C — ► R be a measurable function 
that associates to a function uj' £ C its local time at in the interval (0, 1] 
if u/ has continuous local time. An alternative representation of the process 
in (4.2) is then 

(4.3) Z n = j Lo</> n (u/)Mi(<L/), n>0, 

Jc 

where Mi is an SaS random measure on C with control measure m, and 
(j) : C — ► C is given by 4>{oj') = u/(- + 1). The stationarity of the increments 
of the fractional Brownian motion implies that the map (j) preserves the 
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measure m. Note that (4.3) is a representation of type (2.2) (with both 
a n = 1 and the Radon-Nikodym derivative equal to 1). A conclusion is that 
the flow (4> n ) and the underlying measure space on which (<j) n ) acts are the 
same, independently of the value of a. Therefore, it is sufficient to prove the 
theorem in the case a = 1 , which we will assume until the end of the proof. 
We continue working with the representation (4.2). Note that, by (2.7), 

in 

]T[Z(x,n + l)(u/)-ZOr,n)(a/)] 

n=0 

(4.4) 

= l(x, m + l)(u> ) — > oo as m — > oo 

outside a subset of £1' x K of measure 0. By Corollary 4.2 of [29], this implies 
that the FBM-i^-local time fractional SaS noise is generated by a conserva- 
tive flow. It also, evidently, shows that the kernel in the representation (4.2) 
has a full support. 

In order to prove that the FBM-i/-local time fractional SaS noise is 
generated by a null flow, we will apply Corollary 2.2 of [33] to the obvious 
two-sided extension of the process to (Z n ,n G Z). For reasons of symmetry, 
it is enough to exhibit a nonincreasing nonnegative sequence w n such that 

oo 

(4.5) ^2w n = oo 

n=0 

and 

oo 

(4.6) ^2w n [l(x,n + l)(uj')-l(x,n)(uj')} < oo 

n=0 

for P' x Leb-almost every (u/,x). 

Let w n = (1 + n)~ e with some 1 — H < 6 < 1. Since 6 < 1, the condi- 
tion (4.5) is satisfied. To check (4.6), it is clearly enough to find a strictly 
positive measurable function g such that 

„ oo 

(4.7) E' / g(x) ^ w n [l(x,n + 1){lo') -l(x,n)(uj')]dx <oo. 
Note that 

oo oo - n+1 



E' / 5 (x) V Wn [/( 2 ;,n+l)(u; / )-K^^)(^)]^ = E u; ™ / E 'ff(^(*)) 

Choose g(x) = exp(— x 2 /2) so that for all t > 

1 



K'g{B H {t)) 



f 1 + t 2H a 2)t/2' 



10 



S. COHEN AND G. SAMORODNITSKY 



where a 2 = Var.B#(l). The left-hand side of (4.7) is then 

00 f n+1 dt 00 1 

(48) E Wn 'n (l + W)V2 ^g o W " ( 1+re2 ^ 2 )l/ 2 <00 

by the choice of 6. Hence, (4.6) follows. □ 

Remark 4.2. It follows from Theorem 4.1 that (for 1 < a < 2) the 
FBM-i7-local time fractional SaS motion is different from the linear frac- 
tional SaS motion (or, more generally, from the self-similar mixed average 
processes of [26] ) since the increments of the latter are generated by dissipa- 
tive flows, and it is also different from the real harmonizable fractional SaS 
motion whose increments are generated by positive flows. 

In particular, the FBM-ff-local time fractional SaS noise can be viewed 
as a long memory process; its memory is longer than that of the linear 
fractional SaS noise, but shorter than that of the harmonizable fractional 
SaS noise. Implications of this will be seen, in particular, when we discuss 
smoothness of the sample paths in the next section. This is a reminder that 
very little is determined merely by the exponent of self-similarity for a-stable 
SSSI processes. 

Remark 4.3. It follows immediately from Theorem 4.1 and Theorem 
3.1 in [33] that the FBM-ff-local time fractional SaS noise is (unlike the 
harmonizable fractional SaS noise) ergodic. It is easy to show that it is also 
a mixing process. Indeed, it suffices to show that, for any < a < b and 
e>0, 

lim (P' x Leb){(u/, x) : a < l(x, l)(u/) < b, 

n— >oo 

(4.9) 

l(x,n + l)(u} ) -l{x,n)(J) > e} = 0; 

see, for example, [11] or [30]. Since the left-hand side of (3.3) is finite, we 
see that 

(P' x Leb){(u/, x) : a < l(x, l)(u/) < 6} < 00, 
and so, given 5 > 0, for K large enough, 

(P' xheb)\(uj',x):a<l(x,l)(u}')<b, sup \B H (t)\ > k] < 5. 

I 0<i<l J 

For such K, 

(P' x Leb){(uj',x):a< l(x,l){u') < b,l(x,n + -l(x,n)(uj') > e} 

< 5+ (P ; x Leh){(u',x):\x\ < K,l(x,n + - Z(x,ra)(u/) > e} 

< S + 2KP'(B H (t) e [-K, K] for some n < t < n + 1). 
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Since the last probability clearly goes to zero as n — > 00, we conclude that 
limsupP' x Leb{(u/, x) : a < l(x, < b, 

n— >oo 

l(x, n + l)(u/) — Z(cc, n)(u/) > e} < 5, 
which proves (4.9), since 5 > is arbitrary. 

We close this section by addressing the point mentioned in Remark 3.2. 
Since, in the case a = 1, the exponent of self-similarity of an FBM-ff-local 
time fractional motion does not depend on H, one may suspect that H does 
not change the law of the process itself (up to, perhaps, a multiplicative 
constant). The following result shows that this is not the case, and so the 
parameter H parameterizes an entire family of different 1-stable SSSI pro- 
cesses that does not have either a Levy 1-stable motion or a straight line 
process as a component (indeed, the former would have introduced a dissi- 
pative component to the flow generating the increment process, while the 
latter would have introduced a positive component to that flow): 

Proposition 4.4. Let a = 1 and < Hi, H 2 < 1 with Hi^H 2 . Then, 
there is no constant C such that 

(Y Hl (t),t>0)±(CY H2 (t),t>0), 

where (YHi(t),t > 0) is an FBM-Hi-local time fractional motion with a = l, 
i = l,2. 

Proof. Assume that H\ < Hi- If some C with the above property ex- 
isted, then we could use the fact that the kernel in the representation (4.2) 
has full support and Theorem 1.1 of [29] to connect the kernels with different 
H. Specifically, there would exist measurable maps 

A-.n' xIr|\{0}, 
$i:(1'xIh1, 

such that 

l Hl (x, n + l)(u/) - l Hl (x, n)(J) 

(4.10) = A{u', x)(l H2 ($x(cy , x),n + l)($ 2 (w', x)) 

-/ ff2 ($i(w',x),n)($ 2 (w / ,x))), neN, 

for P' x x Leb-almost every uj' £ f2', x £ R, where we have added subscripts to 
the local times with the obvious meaning, and P^ is the probability measure 
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on f2' corresponding to the fractional Brownian motion with exponent Hi. 
Adding up, we obtain 

(4.11) l Hl (x,n)(u/) = A{J ,x)l H2 ($i(u' ,x),n)($ 2 (w' »), ueN, 

for P^ x Leb- almost every w'6S)',i£l. 

By (2.6), Markov inequality and boundedness of the moments of the local 
time, there is a finite K such that, for every x €M,t > and e, 5 > 0, 

(4.12) P' 2 (l H2 (x,t) > et l ~ H2+s ) < Ke~ 2 r 25 , 
and so by Borel-Cantelli lemma 

(4.13) P' 2 {Ih 2 {x, 2 m ) > £ 2 m(1 - H2+5) infinitely often in m) = 
for every x S R. By Fubini's theorem, 

(4.14) (P' 2 x Leb)(G c ) =0, 
where 

(4.15) G=Uu/,z), lim ^#^# = 

Therefore, in the definition (3.1) of the process (Y(t),t > 0) for H = H2, we 
can restrict the integral from O' x R to G only and then, in (4.10) and (4.11), 
we will have 

(4.16) (<f> 1 (u',x),$ 2 (u',x))eG 

for all lu e fi', x 6 R. This means that for x Leb-almost every w £ 0',iG 
R, we have 

(4.17) lim '-frffff-O. 

Therefore, there is x G R, such that (4.17) holds P^-a.s. 
However, by (2.6), 

P[(l Hl (x, 2 m ) > 2"^-^)) = P ; (l Hl l) > jpC*i-*>-W>) . 

If (5 < -f^2 — -ffi this then gives us 

liminf P' 1 (Z Hl (x,2 m ) > 2 m ( 1 "^ 2 + 5 )) > Pi(M°> 1) > 0) = 1 > 0, 

contradicting (4.17). Therefore (4.10) is impossible and the proposition is 
proved. □ 
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5. Holder continuity. The fact that the local times of the fractional 
Brownian motion are continuous and monotone in the time variable already 
implies that a FBM-i^-local time fractional symmetric a-stable motion with 
< a < 1 is sample continuous (see, e.g., Theorem 10.4.2 of [34]) and the 
same is true for 1 < a < 2 by the mere fact that H' > 1/a (see Theorem 
12.4.1 of [34]). Our goal in this section is to prove Holder continuity of an 
FBM-i/-local time fractional SaS motion for all < a < 2. 



Theorem 5.1. Let (Y(t),t> 0) be an FBM-H -local time fractional SaS 
motion, < a < 2. Then, it has a version with continuous sample paths 
satisfying 

/ \ \Y(t)-Y(s)\ 

(5.1) sup — 73— NN „ , -i ln < oo a.s. 

Remark 5.2. It is instructive to express the Holder continuity state- 
ment in (5.1) in terms of the exponent of self-similarity H' of the FBM-H - 
local time fractional SaS motion and a, which can be done for For 
such a, (5.1) means that an FBM-if-local time fractional SaS motion is 
<i-H61der continuous with any 

, H'-l/a 

5.2 d< . 

1 — 1/a 

Let, for example, 1 < a < 2. Recall that a linear fractional SaS motion 
with exponent of self-similarity H' > 1/a is d-H61der continuous with any 
d < H' — 1/a [35], while a harmonizable fractional SaS motion is d-H61der 
continuous with any d < H' [16]. In particular, an FBM-i^-local time frac- 
tional SaS motion has smoother sample paths than a linear fractional SaS 
motion with the same exponent of self-similarity, and less smooth sample 
paths than a harmonizable fractional SaS motion with the same exponent 
of self-similarity. This is not surprising if one recalls that the increments of 
an FBM-if-local time fractional SaS motion have "stronger dependence" 
than those of a linear fractional SaS motion, but not as strong as those of 
a harmonizable fractional SaS motion. 

Of course, since Theorem 5.1 only provides a lower bound on how smooth 
the sample functions are, the above discussion should be taken with "a grain 
of salt." We conjecture, however, that the upper bound on the Holder expo- 
nent of the FBM-ff-local time fractional SaS motion cannot be improved. 
In the case H = 1/2, this is shown in Remark 5.3 below. 



Proof of Theorem 5.1. We will use a series representation of the 
stochastic integral (3.1) defining the FBM- //-local time fractional SaS mo- 
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tion. In distribution, 

oo 

(5.3) Y(t) = C a Y / G 1 TJ l/a e X i /2a l J (X j ,t), t>0, 

i=i 

where C a is a finite positive constant that depends only on a, where (Gj), 
(Tj), (Xj) and (lj) are four independent sequences such that (Gj) and (X,) 
are i.i.d. standard normal random variables, (Tj) are the arrival times of a 
unit rate Poisson process on (0, oo) and (lj) are i.i.d. copies of the local time 
process of a fractional Brownian motion; see Section 3.10 in [34]. 

Assume that the sequence (Gj) is defined on some probability space 
(fil, ^i,Pi), while the rest of the random variables on the right-hand side 
of (5.3) are defined on some other probability space (0,2,^2, P2), so that 
the FBM-ff-local time fractional SaS motion in the left-hand side of (5.3) 
is defined on the product of these two spaces. Let 

lj(x, t) — lj(x, s) 

^ (t-s)i-*(logl/(t 

0<s<t<l/2 



K j — SUp , 1 itx -WW J — 1)2, 



and notice that, for a fixed 102 £ ^2, the process in (5.3) is centered Gaussian 
with the incremental variance 

00 

E 1 (Y(t) - Y(s)) 2 = C^T- 2 / a e X l/ a (lj(Xj,t) - lj(Xj,s)f 
(5-4) <{clf J Tf a e X V^A(t-s) 2 ^{\o g ^ 



M(u 2 )(t-s)^- H ){\o gT ^j 



2H 



for all < s < t < 1/2, where M is a P2-a.s. finite random variable on 
(^2,^2,^2) (the latter statement follows from the fact that E^K'j < 00). 

Applying now classical results on moduli of continuity of Gaussian pro- 
cesses (see, e.g., Theorem 2.1 of [8]) we obtain that, for P2-almost every 

\Y(t)-Y(s)\ 

sup 1 n — ett; 77 zj 1 1 m < 00 > Pi-a.s. 

o< s < t <i/ 2 (t - S y- H (log i/(t - s)) H +v 2 

s,t rational 

By Fubini's theorem, 

\Y(t)-Y(s)\ „ „ 

sup - -— tt- — — . < 00, Pi x P 2 -a.s., 

o< s <t<i/2 (t~s) L H (\ogl/(t - s)) H + L / 2 



which is equivalent to the statement of the theorem. □ 
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Remark 5.3. It is easy to show that, at least for H = 1/2, the result 
of Theorem 5.1 is "almost" sharp in the sense that there does not exist a 
function g : (0, 1/2) — ► (0, oo) with 

J™ ^(logl/i)!^ = 

and, with positive probability, 

/ - -\ \Y(t)-Y(s)\ 

(5.5) sup - — ^ < oo. 

0<s<t<l/2 9\t — s ) 

Indeed, assume that such a function, in fact, exists. By the zero-one law for 
stable processes; see Section 9.5 in [34]. (5.5) would then hold with proba- 
bility 1. It follows (e.g., by [28]), that we must have 

\l(x,t) -l{x,s)\ 

(5.6) sup r < oo a.s. 

0<s<t<l/2 9{t — s ) 

for almost every Choose x for which (5.6) holds and note that by the 

strong Markov property and [12], 



P' 



\l(x, t) — l(x, s)\ 

SUP r = OO 

0<s<i<l/2 9{t ~ s ) 

>P'(mf{u>0:B 1/2 («)= 2 :(<i)p'( sup IH "'l" °" = °° 



1W sup MfirKMI 

4/ Vo<s<t<l/4 9(t 



= P' ^inf {u > : B 1/2 (u) = x} < ^ > 0, 

contradicting the necessity of (5.5) to hold with probability 1. 

6. Expansion into absolutely continuous terms. The sample paths of 
(measurable) //-SSSI processes are almost never absolutely continuous with 
respect to the Lebesgue measure, the only exception being the case H = 1 
with the process being the straight line process Y(t) = tY(l) a.s. for all t 
(see Theorem 3.3 of [37]). Nonetheless, there is a school of thought view- 
ing nature as "producing smooth objects," with the others being more of 
a mathematical abstraction. In particular, smooth modifications of various 
mathematical models are of interest; see, for example, the "physical frac- 
tional Brownian motion" of [13]. In this section we use the chaos expansion 
of the local times of fractional Brownian motions due to [9] to construct 
an expansion of the FBM-ff-local time fractional SaS motion into a series 
of absolutely continuous SaS self-similar processes, all with the same expo- 
nent of self-similarity as the original process. We first introduce the required 
notation. 
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For a > 0, let p a i denote the density of a zero mean normal random 
variable with variance a 2 . H n is the nth Hermite polynomial 

, . (-l) n fx 2 \ d n ( ( x 2 \\ 

with Hq{x) = 1. Let (fi'^^P') be a probability space supporting a Brow- 
nian motion (W(s),s € R), and let I n be the nth Wiener-Ito integral with 
respect to this Brownian motion. We refer the reader to [23] for information 
on these notions. Finally, let Kjj be the kernel defined by 

K H (t, s) = (t- sf- 1 / 2 -{H- 1/2) jf ( r - s) H ^ 2 (l - (jj J dr 

for < s < t and equal to zero for other values of s,t. Note that, in distri- 
bution, 

(6.1) B H {t) = ( Var ^ (1) ) V2 j*K H (t,s)W(ds), t>0, 

where Ch is a finite positive constant depending only on H (see, e.g., [2]). 

Theorem 6.1. Let (Y(t),t> 0) be an FBM-H -local time fractional SaS 
motion. In distribution, 

OO OO p p 

(6.2) Y(t) = Y,W n (t):= V / / h n {x,t)(u>')M(cLj',dx), t>0, 

where, for n = 0, 1, . . . , 

h n (x,t) = h n (x,t){d) 

(6.3) 

rai/i a 2 = Var.B#(l), and M an SaS random measure on the space fl'xl 
lozi/i control measure P' x Leb. i?ac/i process (W n (t),t > 0) is a self-similar 
SaS process with exponent of self- similarity H' = 1 — H + H/a and has a 
modification with absolutely continuous sample paths. Moreover, the series 
in (6.2) converges in probability. 

Proof. We first check that each (W n (t),t > 0) is a well-defined SaS 
process. Note that 

Jn' jr 

= {aC n I l 2 ) a V' [ h%(x,t)dx 
m 
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< 



< 



E' 



t p s 2 H {x/a) ( x/a 



M 



2i 



a/2 



dx 



p s2H (x/a) 



x/a 



I n (K H (s,-r n )ds^ 
[E'(I n (K H (s,-)® n )) 2 } 1/2 d s y dx. 



By (6.1) 

E'(I n (K H (s,-)® n )) 2 = n\\\K H ( S] 
Therefore, 

E' / h%(x,t)dx< [ 



( n! )a/2 



n \\l=n\\\K H {s,-)\\T =n\C n H (s 2H ) n - 
x/a s 



o 



P s 2h(x/<j) 







Pi 



X 



as 



S H 

x/a 



a 

ds ) dx 



ds 
7h 



dx. 



Observe that the function ip(x) = pi(x) 1 ^ 2 \H n (x)\ is continuous and bounded 
on the entire real line. Therefore, 



PL 



X 



as 



11 



x/a 



ds 



a 



l ' 2 ds_ 

) S H 



for some finite positive c independent of x. Therefore, 

(6.4) E' / h%(x,t)dx < oo, 

and so each (W n (t),t > 0) is a well-defined SaS process. 

The next step is to check that each process (W n (t),t > 0) is self-similar, 
with the exponent of self-similarity given by (3.2). We will use two simple 
scaling facts. The first is simply 

(6.5) K H (au,w) = a H ~ 1/2 K H (u,w/a) 

for all a > and all u, w, and the second is a consequence of the self-similarity 
of a Brownian motion: for any n > l,m> 1, any square- integrable symmetric 
functions fi,...,f m an d a > 0, 



(6.6) 



(I n (fi(-a)),i = 1, . . . ,m) = (a~ n / 2 I n (f l ),i = 1, . . . ,m). 



We assume, for simplicity, that a = 1. 

Now let m > 1, < t\ < ■ ■ ■ < t m , 0i,...,0 m eM and a > 0. We have 



E' 



dx 
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n—anH 



E' 



3=1 



u 



a H u H 



xI n (K H (au,-)® n )du 



dx 



^a—anH 



E' 



in „f . I I H 



a n u n 



a 



ce—anH+H—aH 



E' 



x I n {K H (au,-r n )du 
f tj Pu'^{y) f y 



dx 



3=1 



o u 



u 



11 



xI n (K H (au,-r n )du 



a-an/2+H-aH 



dy 



xI n [K H [u, 



da 



dy 



,a+H-aH 



E' 



dy 



= - logE expjiC^ 2 £ d ja H 'W n (at 3 )j , 

where the fifth equality follows from (6.5) and the sixth equality follows 
from (6.6). This proves the claimed self-similarity. 

Because of self-similarity, it is enough to prove absolute continuity on the 
interval [0, 1]. The proof will be done in three different cases; each consists 
of checking the conditions of Theorem 11.7.4 in [34]. 

If < a < 1 , we need to check that 



E' 



dh n 



dt 



(x,t) 



dt) dx < oo. 



This is, however, an immediate consequence of the computation leading to 
(6.4). 

If 1 < a < 2, we need to check 



(6.7) 



E' 



dh„, 



dt 



(x,t) 



a \ l/a 
dx) dt < CO. 
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We have, for t > 
dh n 



E' 



dt 



x,t) 



dx 



dx 



< t~ anH (-E'l n (K H {t, .)® n f) a l 2 

( P1 (x/t H )\H n ( X /t H )\r 



4-aH 



dx 



h ,-(a-l)H 



for some < b n < oo. Since 



t -(a-l)H/a dt 



< OO, 



(6.7) follows. 

Finally, in the case a = 1, the necessary and sufficient conditions for abso- 
lute continuity are less convenient to check. However, a stronger statement, 
that the process is absolutely continuous with a derivative in L p [0, 1] for 
some 1 < p < 2, requires checking that 



E' 



dh n 



dt 



(x,t) 



1/p 

dt) dx < oo, 



which follows in the same way as (6.4). We omit the repetitive details. 

It remains to prove that the sequence (6.2) converges in probability. By 
Proposition 4 of [9], for every x and t 

oo 

l(x,t) = h n (x,t), P'-a.s. 

n=0 

By definition (3.1) of the process (Y(t),t > 0), it is enough to prove that 



(6i 



E' 



l(x,t) - h n( x ,t) 
n=0 



dx — ► asm-> oo. 



We will estimate the expectation in (6.8) in two different ways. First note 
that, for every m, 



E' 



l(x,t) - Y h n (x,t) 

n=0 



2\ a/2' 



E'l(x,t) a +[E'[j2hn(x,t) 

\ \n=0 



<C a [E'l(x,t) a + (E'l(x,t) 2 ) a/2 ] 
<C a (E'l(x,t) 2 ) a/2 , 
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where C a is a finite positive constant depending only on a and allowed to 
change from place to place. The argument used in (3.3) shows that 

(E'l(x,t) 2 ) a/2 dx<oo. 

Therefore, given e > 0, one can choose M 6 (0, oo) such that for all m > 1, 

-M 



E' 



l(x,t) - ^ h n (x,t) 

n=0 



dx < e + 



E' 



Af 



- ^ h n (x,t) 

n=0 



dx. 



(6.9) 

Next, we estimate the expectation in (6.8) in a different way. Note that 
by the orthogonality of (fo n )'s with different n, 



E' 



l(x,t) - ^ h n (x,t) 



n=Q 



in 



2\ a/2 



< [E'[l(x,t)-J2h n (x,t) 



n=0 / , 

oo \ 2\ «/2 



E'I ^ /i n (x,t) 

V \n=m+l ^ 

' oo 

^ ^'(^(x,*)) 2 

\ n=m+l 



a/2 



and, as in the proof of Proposition 4 in [9], we conclude that there exist 
some 5 m t — ► as m — > oo such that, for all 



E' 



l(x,t) - ^ h n (x,t) 

n=0 



< 5. 



m,i ■ 



Substituting this bound into (6.9), we conclude that 



E' 



l(x,t) - ^ h n (x,t) 

n=0 



dx <e + 2M8 mt . 



Letting first m — > oo and then e — > proves (6.8), and so the proof of the 
theorem is now complete. □ 

Remark 6.2. It is clear from the proof of the theorem that the deriva- 
tive of each process (W n (t),t > 0) is in L p [0, 1] for a range of p > 1 in all 
cases, and not only for a = l. We will not pursue this point here, however. 



7. Convergence of the random reward scheme. In this section we es- 
tablish the limit theorem in the random reward scheme discussed in the 
Introduction. We start by setting up the notation. Let (Wu , k E Z,i > 1) 
be an array of i.i.d. symmetric random variables whose distribution satisfies 
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(1.1). Let {Vu ,k > 1, i> 1) be an array of i.i.d. mean zero and unit vari- 
ance integer- valued random variables, independent of (wi , k € Z,i > 1). 

Let Sn = v} + • • • + Vn , n > 0, be the ith random walk, z = 1,2,..., and 
define for j G Z and n > 1 

n 

(7-1) ^(i,n;i) = E 1 (^ ) =i) 

fc=i 

to be the number of times the ith random walk visits the state j by time 

ra, 4 = 1, 2, Define (p(j, t; i) for noninteger values of t > by interpolating 

linearly between ip(j, n; i) and ip(j, n + 1; i) if n < i <n+l [we use </?(j, 0; i) = 
0]. Notice that the total reward earned by the ith user by time t can be 
written as 



oo 

r(i), 

h = — OG 



(of course, this is really a linear interpolation for noninteger t). The limit the- 
orem below shows that, if both the number of users and the time scale grow 
at an arbitrary rate, then the properly normalized total reward converges 
weakly to the FBM-l/2-local time fractional symmetric a-stable motion. 
This is related to the convergence result in [15] (which allows more general 
random walks) where only one user is present. 

Theorem 7.1. For every sequence (pn) of positive integers with b n — > oo 
we have, as n—>oc, 



((2/C Q ) lA V w F(i),t>0) 



(7.2) 

weakly mC([0,oo)), where (Y(t),t > 0) is the FBM-l/2-local time fractional 
symmetric a-stable motion defined in (3.1) (with the local time being that of 
a standard Brownian motion). Here, aw is the tail weight in (1.1) and C a 
is the stable tail constant given by 

(7.3) C a =(^J x~ a smxdx s j 

Proof. By extending the probability space on which random objects 
are defined, if necessary, we can construct a sequence of i.i.d. standard Brow- 
nian motions (B^ l \t),t >0), i = 1,2, ... , with jointly continuous local time 
processes (l^\x, t), t > 0, t £ R), 2 = 1,2,..., such that, for every T > 0, 



(7.4) sup 

x&R,0<t<nT 



o 
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in probability as n — > oo, i = 1,2,...; see Theorem 1 of [5]. Define, for n > 1, 

(^V ) 7 ^Ife^oo W"n / 

Notice that, for t > 0, 



1 n 

^(*) := ,,i a+ D/ 2 a/ a E fi(i) (M) - x n (t) 

(7.6) 



(^ +1 )/ 2 ) 1 /° «, 



(n^ Q+1)/2 ) 1/a 



n oo 

X 

i=l fc=— oo 



We first prove that, for every t > 0, 

(7.7) £ n (i)^0 in probability. 

For notational simplicity, we prove (7.7) for t = l. 

First, it follows from the tail behavior (1.1) that there is a constant b > 
such that 

(7.8) |W»|< 6(1 + 11^1) 

(in the sense of stochastic comparison), where (Ru ,k G Z,z > 1) is an ar- 
ray of i.i.d. standard SaS random variables. Therefore, by the contraction 
inequality (see Section 1.2 of [18]) we conclude that 

P(\E n (l)\>e) 

( 1 

< 2P 



(nbt +1)/2 f /a 



n oo 

*E E 4 l) a+K 



i=l k=— oo 



(7.9) x ^k,b n] i)-bi^^,l^>e/b 



< 2P\ 



1 



± £ e«f^,6 n ;i)-6y^)f^,l))> e /(2^ 



x 

2=1 fc=— OO 
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+ 2P 



(a+l)/2xl/a 



\{nbh 



n oo 



X E E RVW{k,b nV i)-b]!^(^l))>e/{2b) 



k 



i=l k=—oc 

:=pi(n) +j>2(n), 

where (e^\k GZ,i> 1) is an array of i.i.d standard symmetric Rademacher 
random variables. We need to show that 

(7.10) Pj( n ) — > as n — > oo for j = 1,2. 

We estimate P2(^)- Note that 

p 2 (n)/2 

-p(' 



n oo 



E E <p(k,bn;i)-bX*l® -£=,1 



A: 



1/ct 



M 1} > e/(26) , 



i=l fc=— oo 

and so the statement (7.10) with j = 2 will follow once we show that 



(7.11) 



i 71 OO / h 



1 







in probability as n — > oo. The expectation of the expression in the left-hand 
side of (7.11) is 



.(«+l)/2 

u n fc=-oo 



< 777^1772^ E ^,^l)-^wf-^,l 



,(«+l)/2 

°n fc=— oo 



X 1 



¥>(A;A;l)-£ /2 i (1) 



< 1 



+ 



1 



(o+l)/2 



6 



x £ 



> fc=— oo 



a/2 
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\fe=-oo ^ 



> l 



l-a/2' 



,(a+l)/2 

Un fc=-oo 



X 1 



V{k,b n -l)-b l J 2 l^ 



< 1 



+ 



b 



(«+l)/2 



£ ip{Kb n -A)-b]i 2 i {l) 

. \ k=—oo 

x (2A/W(& n ) + 1) 1 ~ Q/2 



a/2 



x 1 sup 



:=P2i(«) +P22{n), 
where 

M«(m)=maxf sup I^^V^sup \B^(t)\). 

\0<k<m 0<t<l J 

For the second inequality above, we have bounded a sum from above by 
the number of nonvanishing terms times the largest nonvanishing term. A 
similar argument will be used in the sequel without further comment. We 
have 

P2i{n)<^ m E{2M^{b n ) + l) 

bn 



<c 



b l J 2 = cb~ a / 2 ^0 



,(a+l)/2 n n 
On 



as n — > oo. 



Furthermore, 
P22{n) < 

where 



a/2 



(E(2M^(b n ) + l)l(A 1 (b n )>l)) 



Aj(n) = sup 
fcez 



y{k,b n ;i)-bH 2 l® 



l-a/2 



1 



Using Lemma 1 of [15] and the fact that the largest value of a Brownian local 
time at time 1 has all moments finite, the first expectation in the right-hand 
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3/2 

side is bounded above by c6 n • Therefore, 

P22 (n) < c T i IW 6r/ 4 (£;M( 1 )(6„) 3 / 2 )( 2 -«)/ 3 (P(A 1 (6 n ) > l)f~^ 



<c- 



^e /4 (e) (2 - Q)/3 (^(Ai(6 n ) > l))CW2)/3 



(q+1)/2 



<c(P(A 1 (6 n )>l)) 



(l-o/2)/3 







as n — > oo, 



by (7.4) (as always, c is a finite positive constant that may change from 
instance to instance). Therefore, (7.11) holds, and so we have established 
(7.10) for j = 2. The proof for j = 1 is similar. Thus, we have obtained (7.7). 

The next step is to show that the finite-dimensional distributions of the 
process (X n (t),t > 0) in (7.5) converge to those of (Y(t),t > 0). For this, it 
is enough to show that, for every k > 1, < t\ < ■ ■ ■ < tj. and 6>i, G M, 



3=1 



n — ► oo. 



We will see that this is true for k = 1 and ti = 1; the general case is only 
notationally different. That is, we will show that 



, n oo 

(7.12) V V 



i=l k= 



k 



Yd) 



as n — > oo. 



By Theorem 8 in Chapter 6 of [24] it is enough to prove that, for every 
A>0, 



(7.13) 



^ k=— oo 



Urn nP[ £ ( ^, 1 ) > A^V" 



/ l{x,t) a dx\~ 
Jm. 



and 



lim lim sup ■ 



(7.14) 



(n&, 
x £ 



E 



. fc=— oo 



:.l 



xl 



E 



fc=— oo 



<« f2 f ,a 



(we have used the symmetry of M^'s to simplify the conditions). 
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We start by checking (7.13). The first step is to prove that, for every 
A>0, 



(7.15) lim limsupnPf V wPl™ ( l) > A(n&^ 2 ) 1/a 

K->oo n-oo V VA / , 



0. 



By using the contraction inequality, the stochastic comparison (7.8) and the 
notation following it, it is enough to prove that, for every A > 0, (7.15) holds 
with each being replaced by Rjp , and with each being replaced 

by . The two statements are similar; we only present the argument in the 
case of stable weights. In that case, the expression corresponding to that in 
the left-hand side of (7.15) is equal to 



nP\ 



E 

<\k\>KVK 



k 



1 



R? > KnbT) 



i/2UA* 



and, for some positive constant c, this bounded from above by 

k 



E 

\k\>K^K, 



E 

\>KVb^ 



<c\- a E[ supZW(x,l)) P 

\k\>K^/K 



2c\~ 



sup |J3 (1) (a)| > -4= 

0<s<l yJO, 



P[ sup \B {1 \s)\ >x)dx. 

\0<s<l 



Since the final expression converges to as K — > oo, we have (7.15). 

Now fix K and A > 0. The usual "largest jump" large deviations approach 
(see, e.g., [22]) and the continuity of the local time give us that, as n — ► oo, 

np( £ ^ZW(-^,1)>A(^)1/^ 

~np( max W^l^ ( -£=, l) > \{nb]/ 2 ) l l a 



\k\<Ky/b. 



(7.16) 



~ n 



\k\<K^K 



E P(^ (1) / (1) (^=,l)>A(n6y 2 )i/« 
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~ n 



p(V«Z«(-^=,l) >X(nb 1 r / 2 ) 1 ^dx 



K 

K 
A 



(n6y 2 )P(^ fe (1) /( 1 )(y,l)>A(n6y 2 ) 1 /«)dy 



a 



K 



(see, e.g., (2.7) in [31]). 

(7.13) now follows from (7.15) and (7.16). 
To show (7.14), note that 



n 



l/2\2/a 



E 



E 



, k=— oo 



X 1 



E 



<n P 



k=—oc 
oo 



E 



fc=— oo 



< e(n6y 2 ) 1 / Q 



dx. 



Using stochastic domination and contraction principle as above allows us 

(i) 

to replace the random variables Wt in the above expression by SaS ran- 
dom variables and by Rademacher random variables and, as before, we only 
consider the former (because they have heavier tails). In that case, we have 



n 



f p £ s?»i«(-i.,i 

/ oo , 



n / P 



v fc = — oo 



1/a 



> 



dx 



<cb~^E J2 

fc=— oo 



- a ' 2 dx, 



from which (7.14) would follow once we check uniform boundedness of the 
ra-dependent coefficient above. However, this follows from 



.xei 



sup(^(x,l)) a ( 2 sup \B (1) (t)\ + l 



0<<<1 



< oo. 



Therefore, we have (7.14) and, thus, convergence of the finite-dimensional 
distributions in (7.2). 
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It remains to prove tightness. Write, for M > 0, 
1 n 

ly^E^M) 
i=i 

, n oo 

= ; CTi^ E E ^lO^^MnVa^^^i) 

(»Ow J i=l fc=-oo 

(7.17) 

-, n oo 

(wOn J i=l fc=-oo 

:=y n (t) + Z n (t), t>0. 
Notice that 

p( sup |F n (t)| > o) < 1- [P(for all |fc| < M^(6 n ), |W fe (1) | < M^ a f fl )r. 
\o<t<i / 

Since for large n, with a changing constant c, 

P(for all |*| < M«(6 n ), |W fc (1) | < Mn 1/a lh /2a ) 

= E[P(\wi 1) \<Mn 1 Hl/ 2a )} 2M(1) ( bn)+1 

>tf[l-cM- a n-V /2 ] 2M(1)(6,0+1 

> £exp{-cM _a ra _1 &- 1/2 (2M (1) (& n ) + 1)}, 

we obtain, using the inequality e _:c > 1 — cc for x > and maximal inequality 
for martingales, 

p( sup |r n (t)|>o) <l-[l-cM- a n-X l ^E(2M^(b n ) + l)] n 
Vo<t<i / 

< 1 - [1 - cM' a rT l \ n -» 1 - exp{-cM- Q }, 

as n — > oo. 

Since the last expression converges to zero as M — > oo, it follows from 
(7.17) that it is enough to prove that, for each fixed M, the process (Z n (t), < 
i < 1) is tight. 

However, for all < s < t < 1, we have 

£(Z n (t) - Z n ,( s )) 2 = - 1 P^^ld^l < Mn 1 ^)] 

n z ' a ~ L bn 

oo 

x £ ^ (<p(fc,6 n t;l) - tp(k,b n s;l)) 2 . 

k=—oo 
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Since, for large x, 

9 

E[(W 1 (1) ) 2 l(|W 1 (1) |<x)]<4 / yP(Wf> >y)dy <cx 2 ~ a , 

Jo 

we see that, for large n, 

oo 

E(Z n (t) - Z n {s)f < cb^/ 2 E Y, (v{k,b n t;l) ~ <p{k,b n s;l)f 

k= — oo 

< c(t - sf/ 2 

as in the proof of Lemma 7 in [15]. We can now appeal to Theorem 12.3 in 
[4] to prove tightness of the family of the processes (Z n (t),0 < t < 1) and, 
hence, complete the proof. □ 

8. Discussion and possible extensions. We briefly mention several issues 
related to the model constructed in this paper. 

It is clear that self-similar SaS processes with stationary increments could 
be constructed using local times of self-similar processes with stationary 
increments other than fractional Brownian motions. Symmetric stable Levy 
motions with index of stability between 1 and 2 are obvious examples. One 
could also consider additive functionals other than local times. 

For the random reward scheme considered in Section 7, it is clear that, 
in order to obtain in limit FBM-ff-local time fractional symmetric a-stable 
motion with H ^ 1/2, one has to introduce sufficiently long memory in the 
sequence of steps of each random walk (V^ , k > 1). One way to do it is to 
take a stationary integer- valued sequence with slowly decaying correlations; 
alternatively, a certain reinforcement mechanism could be used. This is left 
for a future work. 

It is also instructive to note that, in Section 7, one obtains the same limit 
regardless of how fast the number of users grows. However, if one consid- 
ers instead (as is common in the literature) a fluid input system, where the 
random reward is not gained instantaneously, but, instead, obtained over a 
stretch of time, it is likely that different limits would be obtained, depending 
on the number of users. Possible limits there would, probably, include frac- 
tional Brownian motions, FBM-ff-local time fractional symmetric a-stable 
motions and, perhaps, additional limit processes. 
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